All-Resonant Control of Superconducting Resonators 



Frederick W. StrauclQ 
Williams College, Williamstown, MA 01267, USA 
(Dated: August 20, 2012) 

An all-resonant method is proposed to control the quantum state of superconducting resonators. 
This approach uses a tunable artificial atom linearly coupled to resonators, and allows for efficient 
routes to Fock state synthesis, qudit logic operations, and synthesis of NOON states. This resonant 
approach is theoretically analyzed, and found to perform signficantly better than existing proposals 
using the same technology. 
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Achieving complete control of the quantum state of 
light is a primary goal in the field of quantum optics 
and quantum information [1 . The preparation and sub- 
sequent interaction of individual photons for quantum 
communication and computation in the optical domain 
remains a challenging enterprise [2 . By contrast, exci- 
tations of the electromagnetic modes of superconducting 
coplanar waveguide resonators can be readily prepared 
and manipulated using Josephson junction circuits [3|. 
These microwave photons have recently been proposed 
as key ingredients in a superconducting quantum com- 
puter [4 . Finding the fastest and most efficient way to 
control these modes is an outstanding problem. 

Great progress has been made by using a tunable ar- 
tificial atom, the superconducting phase qubit, to excite 
and transfer excitations to and from superconducting res- 
onators |5]. Subsequent experiments have prepared indi- 
vidual Fock states [6 , their superpositions and en- 
tangled states of two such resonators [8]. Despite this 
significant progress, complete control of these resonators 
[9l [To appears to require a Fock-state-selective inter- 
action. Such number-state-dependent interactions were 
first seen by Schuster et al [11] and exploited for Fock- 
state measurement p2] in transmon qubit devices, while 
Fock states have also prepared using sideband transi- 
tions [13 . These latter experiments used dispersive (off- 
resonant) coupling of the qubit to a resonator, while the 
former experiments utilized a carefully chosen sequence 
of qubit operations performed off resonance and resonant 
qubit-resonator swaps. As resonant interactions are often 
faster than their off-resonant counterparts, an important 
question is whether complete control can be achieved us- 
ing resonant interactions alone. 

In this Letter, I present precisely such an all-resonant 
method appropriate for superconducting resonators. 
This method is shown to be applicable to the synthe- 
sis of arbitrary superpositions of Fock states of one and 
two resonators. By performing all of the important steps 
with the qubit on resonance with the resonator, these 
applications are found to be significantly faster, more 
efficient, and of higher fidelity than previous proposals. 
While this method is designed for superconducting exper- 
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FIG. 1: (Color online) Jaynes-Cummings Hamiltonian tran- 
sitions, (a) Energy level diagram, with the four basic transi- 
tions uJn,-\- (red solid arrows, on left), cjn,- (blue solid ar- 
rows, on right), Cc;n,'s, (blue dashed arrows, right to left), 
and oJn,^ (red dashed arrows, left to right), (b) Transitions 
frequencies cJu/Stt as a function of n, for ujr/27r = 6 GHz, 
g/27i — 180 MHz, and A = 0. (c) Transition matrix elements 
{j\crx\k) due to a drive on the qubit. (d) Transition matrix 
elements {j\x\k), with x = a + a^, due to a drive on the res- 
onator. Dark values indicate larger matrix elements, while 
the corresponding transitions are indicated by the symbols. 



iments with existing technology, I will present the basic 
principles using the broadly applicable Jaynes-Cummings 
Hamiltonian, so that any tunable atom can be used to 
quickly control any harmonic oscillator mode to which it 
is coupled. 

Model. The Jaynes-Cummings Hamiltonian is 



l-L/h = UJqCFj^Cr- 



g (acr+ 



(1) 



where Uq and Ur are the transitions frequencies and cr+ = 
|1)(0| and are the creation operators for the qubit and 
resonator, respectively, and g is the coupling between the 
qubit and oscillator. Famously, this can be diagonalized 
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by the eigenstates 

|n, — ) = cos^n|0, n) — sin6>n|l, n — 1) 

|n, +) = sin(9n|0,n) + cos(9n|l,n - 1) (2) 

where tan(2^^) = 2^y^/A, A = cj^ — cj^, and the \q^n) 
represents a state with the qubit q = {) oi 1 and oscihator 
number n. The corresponding eigenvalues are given by 



En,±/h = ncj^ + i (^A ± VA2 + Ang^ 



(3) 



Note that the ground state has Eq = {) and is given by 
|0) = |0,0). The energy level diagram is shown in Fig. 
1(a), along with several transitions, whose frequencies are 
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Roughly speaking, the transitions uJn,± correspond 
to excitations of the resonator, while uj^^^ corre- 
sponds to rotations of the qubit. This is most 
clearly seen in the dispersive regime A ^ ^, where 
ujn,± ^ UJr ± [^V^ - (2^ + 1)//^^] is the Kerr- 
shifted resonator transition and cj^,^ ^ ujq-\-{2n-\-l)g'^ / A 
is the Stark-shifted qubit transition. However, in the res- 
onant regime, when A = 0, the latter transition exhibits 
a significantly stronger dependence on n, as shown in Fig. 
1 (b), while a new transition 0;^,^ becomes possible by 
driving the qubit by a term HdrWe = ^f{t)o-x' This can 
be seen by the matrix elements of ax and x = a -\- ^ as 
shown in Fig. 1(c) and (d). 

Previous transmon experiments have probed only a 
subset of the allowed transitions in Fig. [l] The dis- 
persive ujn^^ transitions were studied in [11], while the 
resonant uji^± transitions were studied in p!4l. Observing 
the resonant diagonal transitions cJn,^, a novel type of 
sideband transition, would be a further test of the Jaynes- 
Cummings model for superconducting circuits, and offers 
a new path toward Fock state synthesis and control. 

Previous theoretical studies have considered sudden 
shifts of the qubit frequency from the dispersive regime 
to the resonant regime to swap excitations from the qubit 
to the resonator |l,n) |0,n-hl), and with qubit transi- 
tions performed in the dispersive regime. Here I consider 
controlling the system when the qubit and resonator are 
resonant, by driving the drive the cj^,^ and uJn,s, transi- 
tions simultaneously. As will be shown below, these tran- 
sitions allow for fast, high-fidelity control of the state. Of 
course, for A = 0, the eigenstates of the qub it-resonator 
system are entangled. To decouple the resonator, one 



could adiabatically shift the qubit frequency to the dis- 
persive regime. However, for many applications, this pro- 
cedure is unnecessary, so I will focus on controlling the 
system in the resonant regime. 

Fock state preparation. The first task to consider is 
the preparation of a Fock state |0,A/'), which by using 
the adiabatic decoupling reduces to the transformation 
|0) 1^, — )• The simplest means to do so would involve 
N steps, each performing a two-level transition by alter- 
nately driving the transitions cj^,^ and cj^+i^k^, whose 
precise sequence depends on whether N is even or odd. 
However, since these transitions are sufficiently separate 
in frequency, these transitions can be driven simultane- 
ously. By choosing an appropriate set of amplitudes, one 
can perform this transformation in a single step. 

Specifically, I consider driving functions of the form 

N 

f{t) = [Anit) COsiUnt) + Bn{t) Sm{Unt)] , (5) 

n=l 

where the slowly- varying envelope functions An{t) and 
Bn{t) and frequencies ujn are chosen to optimize the 
transition. An analytical solution to this problem can 
be found in the rotating wave approximation, by letting 
An{t) = and Bn{t) — 0, where the set of amplitudes 
are given by 



'2N 



for 



1, 



2v'n(7V + l-n) for 1< n < TV, 
and the frequencies are given by 



, UJjsf- 
, UJjsf- 
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■1,N 



for N odd, 
for N even. 



(6) 



(7) 



This choice of frequencies follows a zig-zag pattern up the 
Jaynes-Cummings ladder, while the amplitudes are cho- 
sen so that, in a rotating frame [TS|, the driving Hamilto- 
nian is equivalent to an angular momentum operator for 
a spin-A^/2 system. This has the effect of driving a per- 
fect rotation from |0) — ) in a time T = tt/I^o- This 
analytical solution to a population transfer problem was 
first studied by Cook and Shore [16], and was exploited in 
the phase qudit experiment [17j; applications to perfect 
state transfer in qubit networks have also been studied 

How does this single-step transition (with simultaneous 
driving) compare with a multi-step transition (with se- 
quential driving)? To meaningfully answer this question, 
we must limit the drive amplitudes Q^ri to common value 
of l^max- The sequential, multi-step transitions would 
then take a time of T^uiti = (V^ + 2(A^ — l))7r/l]max ^ 
2A^7r/l]max (the factors of >/2 and 2 are due to the 
transition matrix elements). For the single-step tran- 
sition, Eq. (6) shows that Q^ri < ^o(^ + 1)7 so that 
^single ^ ^7r/r^max, a factor of two better than the multi- 
step transition. Studies of perfect state transfer [20^ show 
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FIG. 2: (a) Error as a function of the pulse time T for the 
Fock state transitions |0) \N,—), for N = 4 and using 
the same parameters as Fig. 1. The circles are for a numer- 
ically optimized two- quadrature pulses with M — 1 Fourier 
component, while the squares are for M = 3 Fourier com- 
ponents, (b) Time-dependent probabilities | ('i;|^(t)) |^ as a 
function of time t for an M = 3 pulse with T = 50 ns, where 
|^;) = |0,0),|1,+),|2,-),|3,+), and |4,-). 



that this is in fact an optimal solution, given constant and 
bounded amplitudes. 

For the system at hand, of course, there are more than 
transitions involved, opening up a signficant source 
of error. Isolating a single transition in a mult level sys- 
tem has inspired a number of control methods [2TH23] . I 
have extended the two-quadrature approach of [22^ to 
the multilevel transition case studied here by numeri- 
cally simulating the time-dependent Schrodinger equa- 
tion ihdt\^(t)) = (H + HdriveWI^W), 1^(0)) = |0) with 

the following envelope functions 



M 



An{t) = XI ~ cos{27rkt/T)] , 

M 

Bn{t) = X6l^)sin(2^H/T), (8) 

k=l 

for various values of the total time T and number of 
Fourier components M. For each such case, I numerically 
optimized the fidelity = |(7V, -|^(T))p to find opti- 
mal Fourier components and driving frequen- 
cies ujn [15]. Figure |2|a) shows the resulting error 1 — 
as a function of T of such an optimization performed 
for A/" = 4 and for M = 1 and 3. The time-dependent 
probabilities for an optimized transition with T = 50 ns 
and M = 3 Fourier components are shown in Fig. [2|b); 
the optimized frequencies and Fourier components can 
be found in the Supplemental Information. For M = 1, 
the T = 50 ns pulse has T > 0.99, while for M = 3 
extraordinarily high fidelity transitions are possible in as 
little as 10 ns. This result highlights the power of this 
all-resonant approach, being significantly faster than the 
120 ns taken by the Rabi-swap sequence of [6 (which 
used ^]o/27r ^ 25 MHz and g/27r ^ 18 MHz). 



This method has several advantages over alternative 
proposals for resonator control. First, there is no need 
to shift the qubit frequency to exchange quanta between 
the qubit and the resonator; the 0Jn,s, transitions accom- 
plish this directly, much like sideband transitions in ion- 
trap systems. This allows for larger qubit-resonator cou- 
plings. Second, this approach can also be pursued for 
large couplings and photon numbers, whereas dispersive 
manipulations are limited to < ncdt = 4A^/^^. While 
the multilevel structure of superconducting qubits will 
also limit the transitions [2TH23] , this effect is minimized 
when g < |cJoi — ^12!, precisely when the two-level qubit 
approximation is appropriate for the Jaynes-Cummings 
Hamiltonian. Finally, the large frequency separation in 
the resonant regime allows for much faster transitions 
than using dispersive transitions (on the qubit or res- 
onator). That is, to avoid non-resonant transitions one 
should require l^max ^ Acj, where Au is the smallest 
frequency separation in the problem. Using the expres- 
sions of Eq. (4), the frequency separations in the reso- 
nant regime are Auj^^k^ ~ gjv}!^ ^ Acj^ ~ gjv?!^ ^ much 
larger than those in the dispersive regime Acj^ ~ 2g^ j A 
and Aa;± ~ 2g^ / A (since ^ <C A). For the case of 
Fig. [2|b), the Fourier components satisfy |al^^|, \hhn \ < 
Acj/27r ^ 237 MHz for pulse times T > 10 ns [15]. By 
using the transitions with the largest frequency separa- 
tion, the control pulse can be performed as quickly as 
possible. 

This time advantage is particularly important in the 
presence of decoherence. Following the analysis of [24] . 
and letting the qubit and resonator have dissipation times 
of Tq and T^, respectively, the fidelity of this transition 
\s ^ g-^/(^T(?)g-ArT/(2T^) both the single-step pro- 
cedure and the procedure of [6 . By reducing the total 
time of the sequence by a factor of two (or more), the 
procedure described above will achieve higher fidelity. 

Qudit operations. The Fock states of the resonator can 
be used as an effective (i- level system, called a qudit jlO] . 
Arbitrary operations on the qudit can be composed by 
rotations 1Zj^k{^)- 



Hj^kiO) = exp 



(9) 



The criterion for such arbitrary control is a connected 
coupling graph of such rotations The coupling 

graph for the Jaynes-Cummings ladder is shown in Fig. 
|3ja), which clearly satisfies this criterion. Here all op- 
erations can be performed using the transitions cj^^k^, 
which are indicated by the dark edges; the intermediate 
states |n, +) have been omitted for simplicity. As there 
are many possibilities for simultaneous transitions, there 
are many opportunities for qudit logic synthesis. 

As a simple example, consider the rotation 7^o,4(^)- 
This is not directly possible, as indicated by the lack of 
an edge between |0) and |4, — ) in the coupling graph. 
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FIG. 3: (Color online) Qudit operations with the Jaynes- 
Cummings ladder, (a) Coupling graph for the states |0) 
and |n, — ), for n = 1 7. The darker edges corre- 
spond to couplings through the transitions cj^,^^, while the 
lighter edges correspond to the transitions (jJ±. The in- 
termediate states |n, +) have been omitted for simplicity, 
(b) Schematic rotation 7^o,4(^), consisting of the sequence 
7^0,3 (7r)7^3, 4(^)7^0, 3 (vr), whose second step is indicated by the 
curved arrows (see text). 



but this rotation can nevertheless be performed by the 
sequence 7^o,3(7r)7^3,4(^)7^o,3(7r), as illustrated in Fig. 
[sj^b). The first and final rotations use the simultane- 
ous control pulses described above, which implement the 
swap |0) ^ |3, — ), while state |2, — ) remains unaffected. 
This sequence can be generalized to implement any ro- 
tation IZj^kifi) in three steps. This represents a signif- 
icant advance over the scheme presented in [TO, which 
would require nineteen steps (with four in the dispersive 
regime) . 

NOON state synthesis. As a final example, I consider 
the generation of entangled states between two supercon- 
ducting resonators, of the form 



I* 



NOON 



V2 



(|iV,0) + |0,iV)). 



(10) 



Methods to synthesize such "NOON" states have been 
proposed in ^ and ^26^, the latter implemented in [8|. 
Each of these methods use a sequence of N qubit rota- 
tions and swaps. I now show how the approach of |8l [26] 
can be simplified by the single-step Fock procedure. 

This method begins by preparing two qubits in the 
entangled state 

|*o) = ^(|1)^|0>b + |0)^|1)b), (11) 

with each qubit coupled to a resonator {A or B) in its 
ground state. By adiabatically moving qubit A (B) onto 
resonance with resonator A (5), the joint state is mapped 
to 

I^i) = ^(|1,+)a|0)b + |0)a|1,+)b). (12) 



For N even, one can simply drive the transition |1, +) 
|7V, — ) on A and B in parallel, and move the qubits off 
resonance, yielding |^noon) in a single step. For N 
odd, one must first drive the transition \'^^~) 
(which leaves |0) unchanged), and then drive the transi- 
tion |1, — ) |7V, — ), yielding |^noon) in two steps. 

By using the all-resonant transitions, this method out- 
performs previous proposals in two other ways. First, 
this method can be performed much faster than the first 
proposal [9 , which used the dispersive number-state- 
dependent transitions. Second, this method uses a two- 
level system to control the Fock states, as opposed to 
using a three-level system to selectively excite the res- 
onators [HI [26] . As shown in [24] , the use of higher ex- 
cited states leads to a lower fidelity. Thus, this approach 
is faster, achieves a higher fidelity, and uses fewer steps 
than previous proposals. 

Conclusion. There are many applications of entan- 
gled resonators for quantum measurement. Bell inequal- 
ity tests, and quantum information, which have been de- 
scribed elsewhere [24]. The all-resonant approach pre- 
sented here provides an attractive route to enabling such 
applications, using existing technology. 

There are two limitations to this approach, to which 
we now turn. To decouple the qubit-resonator state, we 
might use an adiabatic mapping of the uncoupled qubit- 
resonator state |0,n) to the Jaynes-Cummings states 
|n, — ). One means to alleviate this is to use a three- 
level system such that ujr = 0J12 [3 [TOl and apply res- 
onant pulses to perform the mapping |0,n) \^^~) = 
(|l,n) — |2,n — 1)) / a/2 , followed by a control sequence on 
the states |n, — ). An alternative is to use the Kerr effect 
alone, as recently analyzed in [27 , and drive the cj^,- 
transitions, albeit at a slower rate. The second limita- 
tion is the direct excitation of qubit states, which thwarts 
taking full advantange of superconducting resonators, 
namely their potential for longer coherence times. How- 
ever, current experiments have shown that qubits can 
have coherence times approaching 100 /is [28], so that 
the method proposed here should allow for complex con- 
trol of states with photon numbers of ten or more. 

In conclusion, I have proposed an all-resonant scheme 
to perform a range of quantum control protocols for su- 
perconducting resonators. This approach has many ad- 
vantages over previous studies, and extensions may un- 
lock the power of on-chip microwave photons. Finally, 
many of these ideas can be directly applied to other sys- 
tems of light interacting with real and artificial atomic 
systems. 

I gratefully acknowledge discussions with K. Jacobs 
and R. W. Simmonds. This work was supported by the 
NSF under Project No. PHY-1005571. 
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I 

SUPPLEMENTAL MATERIAL 



Introduction 



In this supplement we detail the multi-frequency control pulses used in the text. The eigenstates of the resonant 
Jaynes-Cummings Hamiltonian (with h = 1) 



are |0,0) and 



|n,±) = -^(|0,n)±|l,n-l)). 



where |g,n) = |g)qubit ® \n) oscillator- For H > 0, the eigenstates have the nonzero matrix elements 



(n + 1, +\(Jx\n, -) = -(n + 1, -\<Tx\n, +) 



2' 



while the n = state has the matrix elements 



{l,±|(j,|0,0) =± 



^/2■ 



(13) 
(14) 

(15) 
(16) 



To use the spin rotation method to implement the Fock state preparation |0) — ) described in the text, we 

first define the states 



{K)} 



|0), |1,-),|2, +),■•■ ,|7V,-) for iV odd, 
|0),|1,+),|2,-),--- ,|iV,-) for iV even. 



(17) 



Using this basis and the driving term (Jxf{t), the time-dependent Hamiltonian takes the form 

f Eo hi ••• \ 
hi El ft.2 

m) 



h2 E2 
\ • 







■ • ■■ hN 

hN En J 



(18) 
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where ■Ho|u„) = and 

..^1|±V^/W forn=l, 
^nW 2\±f{t) forn>l. 

The signs for hn{t) alternate depending on whether TV is odd (—,+,—,•••,—) or TV is even (+,—,+,••• , — ). 



Rotating Wave Hamiltonian 



By performing a transformation to an interaction picture with respect to T-Lq , we find the rotating frame Hamiltonian 
U = e'^^^n(t)e-'^^^ - Ho is given by 



( /iie-^^i^ 
h2e'^^^ 



\ 







V 



/live^'^^* J 



(20) 



where cOn = En - ^n-i- 
For the driving field 



AT 



f{t) = y^Qn COs(cJnt), 

the rotating wave approximation takes the form 



(21) 



hn{t) 



1 f V^lli 
4 t 



for n = 1, 
for n > 1, 



(22) 



where we have used cos(co't) = (e*^^ + e *^*)/2, neglected all terms with e^^^^ ^k)t ^ j ^ ^ ^i(ojj+ojk)t ^ g^j^^^ suppressed 
the alternating signs for the matrix elements. The final rotating wave Hamiltonian is 

/ V^^i ••• \ 
a/20i ^2 : 



^2 '•• 
V 



Qn J 



(23) 



Setting 



— ^0 I 



/2N 



for n = 1, 



the matrix elements reduce to 



2^n(A^ + 1 -n) for 1 < n < AT, 



(n|?^RWA|n+l) = ^V^(7V+l-n). 



(24) 



(25) 



Letting j = A/'/2 and m = — j+n, we have ^Jn{N + 1 — n) = ^ (j + m) (j + 1 — m) , so that Hbwa = no(«^+ + '>^-)/2 
Jcc where J^^ is an angular momentum operator for a particle of spin j = N/2. 



7 



Numerical Optimization 

The numerical optimization results reported in Fig. (2) in the text were obtained by a simplex search in Matlab 
(fminsearch) over simulations of the time-dependent Schrodinger equation to find the frequencies uJn and Fourier 
amplitudes ai^\bi^-' for the envelope functions An{t) and Bn(t). Note that the transition frequencies uJn/'2.'K calculated 
from the eigenstates of Ho (with u;/27r = 6 GHz and g/27: = 180 MHz) are given by {6180, 5565.4, 6566.3, 5328.2} MHz 
. The search results for several values of T are reproduced i n the following Tables: 



TABLE I: Numerical optimization of the N = 4 transition with M = 1, with frequencies in MHz. 



Time (ns) 


n 




ai^V27r bl^^/27r \ -T 


10 


1 


6252.7 


209.2 


61.0 


0.22 




2 


5947.2 


87.5 


-169.2 






3 


6661.5 


59.7 


-zl7.o 






4 


5450.1 


286.4 


-159.1 




20 


1 


6226.9 


87.1 


69.9 


0.11 




2 


5540.1 


162.9 


-64.9 






3 


6572.4 


99.9 


26.1 






4 


5291.5 


167.9 


-13.6 




30 


1 


6200.4 


19.1 


-69.6 


0.031 




2 


5557.3 


56.9 


29.8 






3 


6584.1 


59.0 


-37.9 






4 


5299.7 


32.7 


101.4 




40 


1 


6194.7 


14.9 


-50.9 


0.014 




2 


5567.2 


44.9 


-5.97 






3 


6561.8 


46.0 


27.5 






4 


5307.9 


32.6 


57.9 




50 


1 


6195.3 


7.54 


-52.3 


0.005 




2 


5579.1 


20.7 


-67.1 






3 


6564.4 


29.7 


-0.03 






4 


5330.2 


4.83 


-70.3 





TABLE II: Numerical optimization of the N = 4 transition with M = 3, with frequencies in MHz. 



Time (ns) 


n 










6i^V27r 


a^fV27r 


6^fV27r 


l-J" 


10 


1 


6118.7 


207.9 


74.2 


-896.3 


-102.2 


-194.8 


-86.5 


3.0 X 10~^ 




2 


5679.1 


340.1 


111.3 


240.4 


-369.4 


-6.03 


91.9 






3 


6694.1 


365.9 


-202.2 


55.36 


-282.3 


41.7 


100.7 






4 


5403.4 


336.7 


18.8 


-0.08 


-295.7 


148.5 


517.9 




20 


1 


6248.7 


134.4 


101.0 


76.8 


-28.2 


-36.6 


-5.58 


2.9 X 10~^ 




2 


5437.0 


178.8 


11.8 


-38.3 


33.3 


108.9 


23.9 






3 


6666.4 


121.1 


0.43 


84 


-47.1 


73.9 


-27.4 






4 


5209.8 


251.9 


36.2 


-67.2 


36.4 


4.61 


-6.6 




30 


1 


6193.9 


11.0 


14.3 


39.4 


-22.9 


7.92 


-54.8 


2.3 X 10~^ 




2 


5560.6 


63.1 


43.5 


-12.9 


-9.31 


1.90 


27.0 






3 


6591.5 


82.2 


7.62 


11.7 


-3.17 


-21.2 


-16.5 






4 


5319.6 


60.2 


0.54 


0.26 


-65.9 


41.9 


-33.3 




40 


1 


6162.1 


4.66 


13.2 


28.3 


29.9 


0.06 


-0.13 


5.3 X 10~^ 




2 


5554.3 


32.9 


18.1 


12.7 


-41.0 


-1.32 


-2.50 






3 


6550.0 


48.8 


34.4 


-31.5 


18.7 


5.39 


38.6 






4 


5343.0 


0.38 


-47.9 


30.4 


-45.6 


0.67 


-18.8 




50 


1 


6182.8 


5.33 


-22.6 


0.22 


-7.19 


23.6 


20.1 


7.0 X 10-^ 




2 


5560.8 


29.1 


21.9 


25.3 


7.44 


-15.5 


12.2 






3 


6574.8 


23.5 


-18.3 


-9.12 


4.15 


18.0 


8.40 






4 


5333.3 


4.82 


-51.3 


6.27 


-23.0 


6.07 


-19.0 





